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^ , trix elements, and uses antenna functions to subtract the real radiation singularities. Up to now, NLO 

antenna functions could be derived in a systematic manner only for hard quark-antiquark pairs, while the 
gluon-gluon and quark-gluon antenna functions were constructed from their limiting behaviour. In this 
paper, we show that antenna functions for hard quark-gluon pairs can be systematically derived from 
an effective Lagrangian describing heavy neutralino decay. The infrared structure of the colour-ordered 
neutralino decay matrix elements at NLO and NNLO is shown to agree with the structure observed for 
parton radiation off a quark-gluon antenna. 



Abstract 

The computation of exclusive QCD jet observables at higher orders requires a method for the subtrac- 
tion of infrared singular configurations arising from multiple radiation of real partons. One commonly 
used method at next-to-leading order (NLO) is based on the antenna factorization of colour-ordered ma- 



1 Introduction 



Experimental measurements of jet production observables are among the most sensitive tests of the tlioory of 
Quantum Chromodynamics (QCD), and yield very accurate determinations of QCD parameters [1], especially 
of the strong coupling constant ag. At present, the precision of many of these determinations is limited not 
by the quality of the experimental data, but by the error on the theoretical (next-to-leading order, NLO) 
calculations used for the extraction of the QCD parameters. To improve upon this situation, an extension 
of the theoretical calculations to ncxt-to-ncxt-to-lcading order (NNLO) is therefore mandatory. 

In the recent past, many ingredients to NNLO calculations of collider observables have been derived, 
including the universal three-loop QCD splitting functions [2] which governs the evolution of parton distri- 
bution fimctions at NNLO. The massless two-loop 2^2 and 1^3 matrix elements relevant to NNLO 
jet production have been computed [3] using several innovative methods [4], and are now available for all 
processes of phenomenological relevance. The one-loop corrections to 2 — > 3 and 1 — > 4 matrix elements 
have been known for longer and form part of NLO calculations of the respective multi-jet observables [5,6]. 
These NLO matrix elements naturally contribute to NNLO jet observables of lower multiplicity if one of 
the partons involved becomes soft or collinear [7]. In these cases, the infrared singular parts of the matrix 
elements need to be extracted and integrated over the phase space appropriate to the unresolved configura- 
tion to make the infrared pole structure explicit. Methods for the extraction of soft and collinear limits of 
one-loop matrix elements are worked out in detail in the literature [8]. As a final ingredient, contributions 
from the tree level 2 — > 4 and 1^5 processes also contribute to (2 2)- and (1 3)-typc jet observables 
at NNLO. These contain double real radiation singularities corresponding to two partons becoming simul- 
taneously soft and/or collinear [9, 10]. To determine the contribution to NNLO jet observables from these 
configurations, one has to find two-parton subtraction terms which coincide with the full matrix element and 
are still sufficiently simple to be integrated analytically in order to cancel their infrared pole structure with 
the two- loop virtual and the one- loop single- unresolved contributions. Several methods have been proposed 
recently to accomplish this task [11]. Up to now, only one method has been fully worked through for an 
observable of physical interest: using the sector decomposition algorithm [12, 13] to analytically decompose 
both phase space and loop integrals into their Laurent expansion in dimensional regularization, and subse- 
quent numerical computation of the coefficients of this expansion, results were obtained for e+e~ 2j [14] 
and pp H + X [15] at NNLO. In contrast to all other approaches, in the sector decomposition method 
one does not have to integrate the subtraction term analytically. 

In [16], we described the construction of NNLO subtraction terms for e+e~ 2j based on full four- 
parton tree-level and three-parton one-loop matrix elements, which can be integrated analytically over the 
appropriate phase spaces [13] . Subtraction terms derived from full matrix elements can be viewed as antenna 
functions, encapsulating all singular limits due to unresolved partonic emission between two colour-connected 
hard partons [6,17]. In particular, process-independent antenna functions describing arbitrary QCD multi- 
particle processes can be directly related to three-parton matrix elements at NLO (one unresolved parton 
radiating between two colour-connected hard partons) and four-parton matrix elements at NNLO (two un- 
resolved partons radiating between two colour-connected hard partons). 

Up to now, antenna subtraction terms (at NLO) were obtained by construction (i.e. by inspecting all 
limits they had to contain), in part from the full matrix elements, and in part by using supersymmetric 
(SUSY) relations between matrix elements containing fermions and bosons [6]. A systematic procedure to 
derive antenna functions at NLO and beyond is not available up to now: this paper aims to contribute to 
such a formalism by showing that quark-gluon antenna functions can be derived systematically from physical 
matrix elements obtained from an effective Lagrangian. 

The NNLO subtraction terms derived from four-parton matrix elements with a hard quark-antiquark 
pair in [16] were used subsequently [18] to compute the a^Cf. correction to e+e^ 3j at NNLO. To extend 
this calculation to the remaining colour factors, further subtraction terms must be derived. In particular, the 
subtraction terms of [16] are sufficient only for processes where unresolved partons are radiated from hard 
quark-antiquark pairs: they form the quark-antiquark antenna functions at NNLO. Besides quark-antiquark 
antennae, e"^e~ 3j also contains radiation from hard quark-gluon pairs, the quark-gluon antenna function. 
In the spirit of [16], it should be possible to extract these antenna functions from the matrix elements 
appearing in the NLO and NNLO corrections to a physical one-particle decay process yielding a quark-gluon 
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final state at leading order. It is the purpose of this letter to show that such a process can be described by 
an appropriate colour ordering of the decay of a massive neutralino into a massless gluino and a gluon, and 
to derive the resulting quark-gluon antenna subtraction terms at NLO and NNLO. 



2 Effective Lagrangian and Feynman rules 

To obtain the correct quantum numbers for a quark-gluon antenna function, one has to consider the decay 
of an off-shell spin-1/2 particle into an on-shell spin-1/2 particle (massless quark) and an on-shell spin-1 
particle (gluon). Since the final state quark is in the triplet representation of SU{3), while the gluon is in 
the octet representation, this implies that the initial state spin-1/2 off-shell particle should also be in the 
triplet representation. SU (3) gauge invariance does however forbid external off-shell states. 

In the colour-ordered formulation of QCD tree-level amplitudes [19,20], one decouples the colour quantum 
numbers of the partons from their Lorentz and Dirac structure. Using this formulation, one can in particular 
represent a parton in the adjoint representation as superposition of two partons (with identical momenta) 
in the fundamental representation. It is thus possible to construct the colour ordered quark-gluon antenna 
functions from the SU{S) gauge-invariant decay of an off-shell spin-1/2 singlet state (neutralino) into a 
spin-1/2 octet state (massless gluino) and a spin-1 octet state (gluon), as we shall show below. 

This decay process occurs in the minimal supersymmetric standard model (MSSM, [21]), where it is 
mediated through a loop involving supersymmetric particles. For the purpose of this study, it is sufhcient to 
describe this process through an effective Lagrangian, whose parameters are obtained by integrating out the 
virtual particles in the loop. In the context of the electroweak sector of the MSSM, this effective Lagrangian 
was first derived by Haber and Wyler [22], to describe heavy neutralino decay into a light neutralino and a 
photon. Its generalization to neutralino decay into gluino and gluon is straightforward: 

Ant = z^^a'^Vx^;, + (h.c.) , (2.1) 

which couples a gluino {ipg) and a neutralino {ipx) to the QCD field strength tensor F^^^. The coupling rj 
has inverse mass dimension, and the commutator of the 7-matrices is 

^^'^ = ^[7^7''] . 

It should be noted that this process was discussed previously in the literature in [23], where however no 
effective Lagrangian was stated. 

The Feynman rules following from the Lagrangian (|2.1|) are 



x(p) 




'\ g (ki,£i J 



r(ko) 

fsss. g\K^l,^) = -asvr'^'y,. , (2.2) 

k g'(k2,£2,v) 

where the arrow indicates the direction of fermion flow. It should be noted that Majorana particles also have 
a fermion flow direction, which does however not coincide with the fermion number flow. The momenta are 
always incoming. 
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Besides these Feynman rules and the standard QCD Feynman rules, one needs moreover the Feynman 
rule for the gluon-gluino-gluino coupling [24] 



r(ko) 

= -gsF'^r ■ (2.3) 

rcki) 

The effective coupling 77 can be computed in the MSSM, it was discussed in [22,23]. In the present 
context, its value is irrelevant, but we do have to take into account that rj is renormalized at one loop. Its 
QCD renormalization constant reads 

a, 1 //3o , 3A^ 




^.-l-^-^y + -j+C^(«f), (2.4) 
with 

6 

In SUSY QCD, both Z,, and f3o are modified. Throughout this paper, we systematically ignore the SUSY 
QCD corrections and restrict ourselves to the subclass of contributions which preserve the QCD renormal- 
izations l|2.4|l and l|2.5|l . 



3 Colour-ordered amplitudes in neutralino decay 

The basic process for the decay of a neutralino into a gluino plus partons is x{q) giPi)giP3)- Its amplitude 
rScLcls 

z77<5"^"W?^(pi,P3) . 

To display the colour-ordered structure of this amplitude, and to illustrate the relation to quark-gluon 
amplitudes, we multiply it with \/2T°'_^\^ [20]. In the squared amplitude, this factor corresponds to inserting 
unity since 

V2T^\^V2T^X = -^"'"^ • 

The resulting amplitude is then 

-^^3 = ^vV2T:xKiP^^P3) ■ (3-1) 
From this structure, it can be seen that the amplitude contains two colour connected (hard) partons and 
therefore two antennae: 

1. A quark-gluon antenna, with quark momentum pi, quark colour index ii, gluon momentum p^ and 
gluon colour index 03. 

2. An antiquark-gluon antenna, with antiquark momentum pi, antiquark colour index 22, gluon momen- 
tum p3 and gluon colour index 03 . 

This derivation is displayed pictorially in Figure ^ It becomes evident that the Majorana nature of the 
gluino allows it to represent both a quark and an anti-quark. 
The squared matrix element is 

T^.il') ^ \MIJ' = {N' - 1) \M9JP1,P3)\' = 4 (iV^ - l) r^^l ■ (3-2) 

T~'g{q^) serves as normalization for antenna functions obtained from higher order corrections to this matrix 
element. 
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^3 ■S^^'^'T^, = 





Figure 1 : Colour flow contained in tree level decay x 
tal) colour indices. 



gg. Double (single) lines denote adjoint (fundamen- 



To demonstrate the cancellation of infrared divergences at NLO, we compute the renormalized one-loop 
QCD correction to the x{q) — > 9{Pi)9{Ps.) decay, 
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12 



+ 



-I+3C3). 



737r4\ 
1440y 



(3.3) 



The infrared poles of this one-loop correction can be expressed in terms of the infrared singularity operator 
[25] 



N 



1 



3 /3o 
4e 2Ne 



as 



2r(l - e) 

Voles {Ti^iq^))^{l^) 4Relg\e,q^)TS>^iq^) 
r(i) 



(3.4) 



(3.5) 



This expression has to be compared to the 2ReJ^^ (e, q^), which is obtained in the decay of a virtual photon 
into a quark-antiquark pair 7* qq at one loop [16]. The factor 4 in H3.5|) appears since the leading 
order process X ~^ 99 contains two distinct quark-gluon antennae, in contrast to the single quark-antiquark 
antenna contained in 7* — > qq. 



4 NLO antenna functions 

Two different emissions off a quark-gluon pair appear at NLO: either the emission of an additional gluon 
or the splitting of the gluon into a quark-antiquark pair. In the context of the ncutralino decay, these 
correspond to the tree level processes x ~^ 999 f^nd x ~^ 599- 

The tree level amplitude for xio) 9{Pi)9{P3)9(P4:) contains only a single colour structure, y"i"3a4 j-^-^ 
order to relate this colour structure to the colour-ordered quark-gluon antennae, we multiply with ^/2T°'\ : 



(4.1) 



showing that the two colour-ordered amplitudes in this matrix element (corresponding to the two different 
orderings of the gluons along the quark-antiquark iii2 colour line) are equivalent to each other up to an 
overall sign because of the identical momenta of quark and antiquark. Squaring the matrix element and 
dividing by a symmetry factor to account for identical gluons in the final state yields 



l-^LasJ' = V'a' {N' -l)N -|A4^(pi,P3,P4)| 



(4.2) 
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with 



M?^^(pi,P3,P4)|' = 4(l-e) 



2S134S14 , 2sfg^Si3 (l-e)si34S34 (l-e)Si34S34 



S13S34 



S14S34 



Sl3 



+ 



Sl4 



_^2£i3£i4 _^ g^^^^ ^ (1 _ ^ 5^^) 

S34 



8s 



134 



(4.3) 



The behaviour of this matrix element in the kinematical hmits where one parton becomes unresolved is 
as follows: 



1. CoUinear limits: 



|2 91II93 



i2 Bl\\9i 



|2 S3 1134 



Sl3 
§14 

(47ra,) 7:0(si34)-^iVP„(z), 



(4.4) 



with ^; being the momentum fraction of one of the collinear partons and the splitting functions 



1^19 



1 + 
1-z 



- e(l - z) , 



Pn. 



l-z 



+ z{l-z) 



2. Soft limits: 



2513 

S14S34 



S1S3S4 I 

|2 S4->p 



S3S4I 



(47ra,) 7^';(si34)iV 
(47ra,) 7^°(si34)A^ 



(4.5) 



Comparing these limits to the limits of colour-ordered QCD matrix elements, one observes that the collinear 
q ^ qg limit contains a colour factor N/2 in QCD, while the collinear g ^ gg limit derived here contains 
a colour factor N. This is precisely what was expected from the discussion in the previous section, since 
the neutralino decay matrix element considered here contains both a quark-gluon and an antiquark-gluon 
antenna. On the other hand, the collinear g ^ gg limit appears here with the same colour factor as in colour 
ordered QCD matrix elements, indicating that the collinear g gg limit is to be split between both antenna 
functions, as discussed in [6,17]. Finally, the matrix element derived here contains two soft limits with the 
soft eikonal factors as expected in QCD, again reflecting the presence of two antennae. 
Integration over the dipole phase space [13] yields 



d$ 



1 



D,ggg 2l"'"9i93S4i 



2 
72 



421 1197r2 
~^ 18~ 



10 3 4 77r2 
3^ ^ y ~ ~6~ 

250 _ TItt^ 
~9~^^ ~ 720 



209 357r2 
~6 



50 



18 



Ca U 



(4.6) 



The tree level amplitude for x{q) — » 5(^1)9(^3)^(^4) contains only a single colour structure T""^]^, which 
is again contracted with V2Ti"i2 



Ml 



(4.7) 
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yielding 

\^l,s,f = V'9'^^\MI,{PUP3,P,)\' , (4.8) 

with 

\MlM^p„p,r = 4 (1 - .) f 2^^ii±^ + 2 (.13 + .14) ] - 16^ . (4.9) 

The only singular configuration contained in this matrix element is the collinear quark- antiquark limit, which 
is as follows: 

\^Ls,f (4'^«^) r^^isu,)—P,^M , (4.10) 

S34 

with the collinear splitting function 

2z(l^z) 

Integration over the dipole phase space [13] and summing over final state quark flavours yields 



^-1+ -3 + - -9+- + -C3).^ + 0(.^) 



(4.11) 

Summing over both three parton final states, we find 

Voles {T^^^iq') + T^^,iq')) = - (g) 4Re/W(6, q') T§^{q^) , (4.12) 

such that the NLO corrected neutralino decay rate into gluino plus partons is finite: 

Voles {T^^iq')) + Voles {T^,g{q^) + T^,,{q^)) - . (4.13) 

It has to be emphasised in this context that we considered here only the QCD corrections to the neutralino 
decay, not the SUSY QCD corrections. At NLO, inclusion of SUSY QCD corrections would both modify 
the renormalization H2.4|) of the effective coupling 77, and include a real radiation contribution from the three 
gluino final state (which has however the same singularity structure as the gqq final state). The omission 
of these corrections is deliberate, since we want to derive the QCD quark-gluon antenna functions. In the 
following section, we demonstrate that the NNLO infrared singularity structure of QCD quark-gluon antenna 
functions is reproduced correctly by the QCD corrections to neutralino decay into gluino plus partons. 

5 Structure of NNLO antenna functions 

In the NNLO calculation of jet observables, two different types of antenna functions are required: (a) the 
one-loop correction to the three-parton antenna functions which appeared at NLO in tree-level form, and (b) 
the tree- level four-parton antenna functions. In this section, we present all neutralino decay matrix elements 
needed for the derivation of these antenna functions, and demonstrate that these matrix elements contain 
the same infrared singularities as processes involving final state emission off a quark-gluon antenna. 

The renormalized one-loop corrections to the three-parton antenna functions have the same colour struc- 
ture as their tree level counterparts listed above. In their computation, closed gluino loops are omitted, since 
these form part of the SUSY QCD corrections. Consequently, renormalization of the coupling constant is 
done using the QCD /3-function. To expose the infrared structure of the resulting one-loop matrix elements, 
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they are integrated over the corresponding dipole phase space [13], yielding 

1 



+- 1 — 



+NNf\ 



20977 2097r2 



108 



12 



9 

2? 



56 
" 3? 

19499 41957r2 



1 / 1835 7l7r2 



36 



27 



72 



12 



695 
— C3 



9957r'' 
720 



4 20 


1 /275 


77r2\ 


/287 


357r2 


100 


3e3 ^ 9e2 


^7 36 - 


9 j 


Hi2- 


27 


9 



J d<f>c,3,,- 2Re (Xg,,3,-,A1 



1,* 

ff 193(74 



3e3 



67 
18^ 



1 /326 
e 



+ 



1 


35 


1 


6^ 


36e2 ^ 


e 


1 


/91 


7r2 


9?" 




27 



509 7r2 
108 ^ T 



27 9 
1670 357r2 



Oie) 



'9215 2757r2 94 
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(5.1) 



V 216 



72 



<3 



81 



31 

2r + y^^ 



(5.2) 



Two different four-parton final states appear in the quark-gluon antenna functions at NNLO: qggg and 
qq'q'g. The Lorentz and Dirac structure of these antenna functions is contained in the neutralino decay 
processes x ~^ 9999 ^^nd x gqqg- In contrast to the tree level three-parton neutralino decay matrix 
elements, which contained only one non-trivial colour ordering each, these four-parton matrix elements both 
contain several colour-orderings. 

To expose the colour-ordered subamplitudes contributing to x{q) ~^ 9{pi)9{P3)9{Pi)9{P5)i we again 
contract the amplitude with \/2T°'\ . The amplitude can then be expressed as sum over the permutations 
of the gluon colour indices: 



M 



91939^95 



iTjg 



1 

V2 



E 



(ij\fe)ePc (3,4,5) 



1 

N 



M~gggipi, P^,PJ,Pk) , (5.3) 



where the sum runs only over cyclic permutations, since the colour-ordered amplitudes ^l^^ggg each contain 
the difference of two colour-orderings which are inverse to each other, as shown in Figurc|2] It can be shown 



that the 1/A^-term in the above expression does not contribute to the physical scattering amplitude [20]. 
The resulting squared matrix element, averaged over identical final state gluon permutations is 



3! 



\M°~ I 
r 51S3S4S5 I 



, , 7V2 ~ 1 1 „ 

772 q4 

' ^ 16 3! 



E 

(i,j,fc)ePc (3,4,5) 



M~gggipl,Pt,Pj,Pk)\ 



(5.4) 



It should be noted that this squared matrix element contains only the leading colour term obtained from 
the squares of the individual colour-ordered amplitudes, as expected in the colour ordered formulation for a 
process with three gluons [19,20]. 

The tree level amplitude for x{q) ^ g{Pi)q{P3)q{Pi)g{Pb) i contracted with ^/2T°'_^^^ contains three colour 
structures. 



M°- - 

91939435 



^i3«2 M^gqqg {Pl , P3, P4, Pb) + '^^i ^4 ^^999 (Pl ' P4 , P3 , P5 ) 

~ J} ^iX M~qqg (pi , P4 , P3 , Ps) 



(5.5) 



7 



il O-i 0-j O-k 12 i2 '^i % ii 



Figure 2: Colour flow contained in the colour ordered amplitude M^ggg(pi,pi,pj,pk) contributing to the tree 
level decay x ^ 5555- 

il 0-5 ii is «2 «3 "^5 ii il ^2 



Figure 3: Colour flow contained in the colour ordered amplitudes M^g^g{pi,p3,p4,p5) (left) and 
M^qgg{pi,P3,P4:,P5) (right) contributing to the tree level decay x 3995- 



The relation between leading and subleading colour ordered amplitudes is 

^Sqqg (Pl : P3 , P4 , P5 ) + M^qqg {Pl,Pi,Pl.,Pb) = ^^qgg {Pl, P3, PA, Pb) 

The squared matrix element reads 



(5.6) 



\^lq^<i4gf = ^'5' {N^ - 1) Nf\ N |A4,-g(pi,P3,P4,P5)|' + \Mlqqg{Pl,Pi,P3,P5)\ 



1 

"n 



M~qqg{pi,P3,P4,P5) 



(5.7) 



It can be seen that this neutralino decay matrix element contains the same colour-ordered antenna 
structures, displayed in Figure 13 as the five-parton matrix element 7* qqq'q' g [9], relevant to e+e" 3j 
at NNLO: gluon (ps) emission between the colour-connected pairs {pi.Ps) and [pi.Pi) at leading colour, and 
gluon emission inside the (p3,P4) pair at subleading colour. In the latter case, the {p3,P5,Pa) system forms a 
colour singlet, such that the gluon p^ acts as a photon and pi becomes a photino which decouples completely 
from any singular limit. 

The four-parton tree-level neutralino matrix elements can be integrated over the tripole phase space [13], 
thus making their infrared singularity structure explicit. 



9999 



d$ 



T,gggg 



3! I »^ 
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4^ 
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1 / 28319 557r2 



e V 144 



188 \ / 2201527 5297r2 
~3~^V ^ V 2592 8~ 



722 SllTT-* 

C3 H 

3 ^ 720 



T?- (q^) = 
.999.9 ~ 



d$ 



T,gqqg 
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\M°~ - I 



n9ii')iir'' 



NNf\ 



5 17 

6e3 ^ 4? 



2239 57r2 
"TO8" ^ ~ 




-OU) 



,(5.8) 
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Nf I 35 1 / 1045 
'"/V" I 6e3 ^ 36^ 7 1^16" ^ T 



/ 28637 357r2 40 
V 1296 ^ 24 ^ 



(5.9) 



The sum of all NNLO subtraction terms yields the following infrared pole structure, which can be 
expressed in terms of NNLO infrared singularity operators [25], 



113 
12^ 



61 
36^ 



55 57r2 
54 ^ T2" 



1 / 27 97r2 



1 /1049 llTT^ 28 



N \ 8e 



2tt 



Re 



?9 



e V 432 



./3o 



3 



+ yC3 



+4lW(e,9')^U9') + 2 



99 



99^ 



(5.10) 



(5.11) 



where (3o is the first term of the QCD /3-function (|2.5(l and the constant K 



(5.12) 



This structure coincides precisely with the singularity structure predicted in [25] for the purely virtual (two- 
loop times tree plus one-loop self-interference) NNLO corrections to a tree level process containing two 
quark-gluon antenna functions. The final state dependent constant Jf^'^' (e, g^) contributes only at 0{e^^): 



J^nci (e,g ) 



4er(l - e) 



(5.13) 



It can be related to the known constants determining the e ^ poles of four-parton two-loop matrix elements 
involving i quarks and j gluons [3] : 



4er(l-e) 



(5.14) 



with 
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H 



(2) 



1. 5 

-Cs H 1 

2^-^ 12 144 



+ —Nl 
27 ^ 



1 



7V2 



1. 41 tt" , 

-C3 + h — A^-" 

4^"^ 108 96 



-C3 + — + -) (2A^^)+f^ 
4^-^ 108 96 / ^ ' ^ 



^ 89 \ TVf 

iVA^F -, 

72 108 / 4iV ' 



3, 3 

-C3 H 

2^"^ 32 



25 \ 2(7V2 _ i-jj^^ 
216 / iV 



48 ~ 216 
1 



NNf 



n-y + iC3) (2A^^) . (5.15) 



In these equations, we decomposed the H-'^^ according to their colour structures. The coefficient {N^ + 1) in 

(2) 

front of the subleading colour contribution to Hq arises due to the fact that abelian diagrams contributing 
to qgggq final states carry this colour structure, such that the generic (planar) leading colour contribution 
is given by just the first term in Hq (see also Equation (3.6) of [9]). 
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It can be seen that contains twice the leading colour and the flavour dependent terms of Hq . The 

(2) 

subleading colour term is absent, and the last term can be identified with the contribution to Hq from 
singularities arising from final states containing two quark-antiquark pairs of identical flavour (Equation 

(4.51) of [16])i_ 

Equations H5.11|l and (|5.15() demonstrate that the NNLO three and four parton contributions to neutralino 
decay into a gluino and massless partons display the same singularity structure as final state observables 
containing adjacent quark-gluon pairs, provided that the colour factors are adjusted correctly. It is there- 
fore possible to construct colour-ordered quark-gluon antenna functions from the neutralino decay matrix 
elements derived here from the effective Lagrangian density ()2.1|) . 

6 Conclusions and Outlook 

QCD antenna functions describe the behaviour of QCD matrix elements in their infrared singular limits, 
corresponding to soft or coUinear parton emission. They are constructed so that they describe all singular 
limits arising from emission of unresolved partons in between the two colour-connected hard partons that 
define the antenna. The quark-antiquark antenna function is directly related to the physical matrix elements 
for 7* — > gg-|-giuons. However, up to now, the NLO quark-gluon and gluon-gluon antenna functions [6] were 
constructed by starting from the quark-antiquark antenna function and adding terms to match the remaining 
limits contained in them. It does not appear feasible to extend this procedure to higher orders. 

In this paper, we demonstrated that quark-gluon QCD antenna functions to all orders can be derived 
from an effective Lagrangian describing the decay of a massive neutralino into a massless gluino and the 
gluon field. In the colour ordered formalism underlying the antenna functions, the Majorana nature of the 
gluino allows it to appear simultaneously as quark and as antiquark. We demonstrated that the physical 
neutralino decay matrix elements reproduce the singular structure of QCD quark-gluon antenna functions 
at NLO and NNLO. We extracted the infrared structure for decay kinematics, as required for jet observables 
without partons in the initial state. By analytic continuation, the matrix elements derived here can also 
be continued to production (leading order process contains partons only in the initial state) or scattering 
(leading order process contains partons in initial and final state) kinematics, where they have to be integrated 
over the appropriate phase spaces. 

All QCD antenna functions can be derived (as opposed to constructed) from physical matrix elements: 
quark-antiquark antennae from the decay of a virtual photon into partons, quark-gluon antennae from 
neutralino decay into gluino plus partons and finally giuon-gluon antennae [26] from Higgs boson decay into 
partons through the effective Lagrangian [27] coupling the Higgs field to the gluonic field strength tensor. 
The NNLO antenna subtraction functions obtained through this procedure will be reported in a subsequent 
publication [28]. 

Acknowledgements 

The authors would like to thank Werner Porod and Daniel Wyler for numerous discussions on neutralino 
interactions. 

This research was supported in part by the Swiss National Science Foundation (SNF) under contracts 
PMPD2-106101 and 200021-101874, by the UK Particle Physics and Astronomy Research Council and by 
the EU Fifth Framework Programme 'Improving Human Potential', Research Training Network 'Particle 
Physics Phenomenology at High Energy Colliders', contract HPRN-CT-2000-00149. 

References 

[1] R.K. Ellis, W.J. Stirling and B.R. Webber, QCD and Collider Physics, Cambridge University Press 
(Cambridge, 1996); 

^In neutralino decay, this identical-flavour contribution corresponds to final states with three gluinos and one gluon. Being 

(2) 

a SUSY QCD correction, it is discarded here, and has to be subtracted from H- ' . 



10 



G. Dissertori, I.G. Knowles and M. Schmelling, Quantum Chromodynamics: High Energy Experiments 
and Theory, Oxford University Press (Oxford, 2003). 



[2] S. Moch, J.A.M. Vermaseren and A. Vogt, Nucl. Phys. B 688 (2004) 101 
A. Vogt, S. Moch and J.A.M. Vermaseren, Nucl. Phys. B 691 (2004) 129 



hep-ph/0403192 



hep-ph/0404111 . 



[3] Z. Bern, L.J. Dixon and A. Ghinculov, Phys. Rev. D 63 (2001) 053007 hep-ph/OOlOOTS ; 

C. Anastasiou, E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeomans, Nucl. Phys. B 601 (2001) 
318 hep-ph/0010212 ; 601 (2001) 347 hep-ph/0011094 ; 605 (2001) 486 hep-ph/0101304 ; 
E.W.N. Glover, C. Oleari and M.E. Tejeda-Yeomans, Nucl. Phys. 605 (2001) 467 hep-ph/0102201 ; 
C. A nastasiou, E.W.N. Glover and M.E. Tejeda-Yeomans, Nucl. Phys. B 629 (2002) 255 
|hep'-ph/0201274 ; 

E.W.N. Glover and M.E. Tejeda-Yeomans, JHEP 0306 (2003) 033 |hep-ph/0304T69| ; 

E.W.N. Glover, JHEP 0404 (2004) 021 hcp-ph/0401119 ; 

Z. Bern, A. De Freitas and L.J. Dixon, JHEP 0109 (2001) 037 | hep-ph/010 9078|; JHEP 0203 (2002) 
018 hep-ph/0201161 ; JHEP 0306 (2003) 028 hep-ph/0304168 ; 

Z. Bern, A. De Freitas, L.J. Dixon, A. Ghinculov and H.L. Wong, JHEP 0111 (2001) 031 
rhep-ph/0109079 ; 

T. Binoth, E.W.N. Glover, P. Marquard and J.J. van der Bij, JHEP 0205 (2002) 060 hep-ph/0202266 ; 

L.W. Garland, T. Gehrmann, E.W.N. Glover, A. Koukoutsakis and E. Remiddi, Nucl. Phys. B 627 

(2002) 107 hep-ph/0112081 and 642 (2002) 227 hep-ph/0206067 ; 

S. Moch, P. Uwer and S. Weinzierl, Phys. Rev. D 66 (2002) 114001 hep-ph/020 7043| ; 

T. Gehrmann and E. Remiddi, Nucl. Phys. B 640 (2002) 379 hep-ph/0207020| . 

[4] F.V. Tkachov, Phys. Lett. lOOB (1981) 65; 

K.G. Chetyrkin and F.V. Tkachov, Nucl. Phys. B192 (1981) 159; 

V.A. Smirnov, Phys. Lett. B460 (1999) 397 hep-ph/9905323 ; 

J.B. Tausk, Phys. Lett. B469 (1999) 225 hep-ph/9909506 ; 

T. Binoth and G. Heinrich, Nucl. Phys. B 585 (2000) 741 hep-ph/0004013 ; 

S. Laporta, Int. J. Mod. Phys. A 15 (2000) 5087 hep-ph/0102033 ; 

T. Gehrmann and E. Remiddi, Nucl. Phys. B 580 (2000) 485 hep-ph/9912329 ; 

S. Moch, P. Uwer and S. Weinzierl, J. Math. Phys. 43 (2002) 3363 hep-ph/0110083^; 

S. Weinzierl, Comput. Phys. Commun. 145 (2002) 357 math-ph/0201011 ; 

C. Anastasiou and K. Melnikov, Nucl. Phys. B 646 (2002) 220 hep-ph/020700i; 

V.A. Smirnov, Evaluating Feynman Integrals, Springer Tracts of Modern Physics (Heidelberg, 2004). 



[5] L.J. Dixon and A. Signer, Phys. Rev. Lett. 78 (1997) 811 |hep-ph/ 9609460 ; Phys. Rev. D 56 (1997) 

4031 hep-ph/9706285 ; 

Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 79 (1997) 3604 'hep-ph/9707309'; 
S. Weinzierl and D.A. Kosower, Phys. Rev. D 60 (1999) 054028 hep-ph/990127^; 
W.B. Kilgore and W.T. Giele, Phys. Rev. D 55 (1997) 7183 hep-ph/9610433 ; 
Z. Nagy and Z. Trocsanyi, Phys. Rev. Lett. 87 (2001) 082001 hep-ph/0104315 ; 

Z. N agy, Phys. Rev. Lett. 88 (2002) 122003 hep-ph/0110315 ; Phys. Rev. D 68 (2003) 094002 
lhep^ph/0307268 



J. Campbell and R.K. Ellis, Phys. Rev. D 65 (2002) 113007 'hep-ph/0202T76 



[6] J. Campbell, M.A. CuUen and E.W.N. Glover, Eur. Phys. J. C 9 (1999) 245 |hep-ph/9809429 



[7] Z. Bern, L.J. Dixon, D.C. Dunbar and D.A. Kosower, Nucl. Phys. B 425 (1994) 217 |hep-ph/9403226 
D.A. Kosower, Nucl. Phys. B 552 (1999) 319 hep-ph/9901201 ; 
D.A. Kosower and P. Uwer, Nucl. Phys. B 563 (1999) 477 hep-ph/9903515|j 



Z. Bern, V. Del Duca and C.R. Schmidt, Phys. Lett. B 445 (1998) 168 hep-ph/9810409'; 
Z. Bern, V. Del Duca, W.B. Kilgore and C.R. Schmidt, Phys. Rev. D 60 (1999) 116001 hep-ph/9903516 ; 
S. Catani and M. Grazzini, Nucl. Phys. B 591 (2000) 435 hep-ph/0007142 ; " ' 

D.A. Kosower, Phys. Rev. Lett. 91 (2003) 061602 hep-ph/03 01069| ; 
S. Weinzierl, JHEP 0307 (2003) 052 hep-ph/0306248. . 



11 



[8] Z. Bern, L.J. Dixon and D.A. Kosower, Phys. Rev. Lett. 70 (1993) 26 77 |hep -ph/9302280|; 
Z. Kunszt, A. Signer and Z. Trocsanyi, Phys. Lett. B 336 (1994) 529 |hep-p h/9405386 ; 
Z. Bern, L.J. Dixon and D.A. Kosower, Nucl. Phys. B 437 (1995) 259Tli^ph/9409393 ; 

E. W.N. Glover and D.J. MiUer, Phys. Lett. B 396 (1997) 257 hep-ph/9609474 ; ^ 

Z. Bern, L.J. Dixon, D.A. Kosower and S. Weinzierl, Nucl. Phys. B 489 (1997) 3 | hep-ph /9610370|; 
J.M. CampbeU, E.W.N. Glover and D.J. Miller, Phys. Lett. B 409 (1997) 503 hep^pfe79706297 ; 
Z. Bern, L.J. Dixon and D.A. Kosower, Nucl. Phys. B 513 (1998) 3 hep-ph /9708239, . 

[9] A. Gehrmann-De Ridder and E.W.N. Glover, Nucl. Phys. B 517 (1998) 269 'he p-ph/970722^ ; 
J. Campbell and E.W.N. Glover, Nucl. Phys. B 527 (1998) 264 ^ep-ph/9710255r 

[10] S. Catani and M. Grazzini, Phys. Lett. B 446 (1999) 143 |hep-ph/9810389 ; Nucl. Phys. B 570 (2000) 
287 hep-ph/9908523 ; 

F. A. Berends and W.T. Giele, Nucl. Phys. B 313 (1989) 595; 

V. Del Duca, A. Frizzo and F. Maltoni, Nucl. Phys. B 568 (2000) 211 |hep-ph/9909464 . 



[11] D.A. Kosower, Phys. Rev. D 67 (200 3) 116003 |hep-ph/0212097| ; 
S. Weinzierl, JHEP 0303 (2003) 062 |hep-ph/0302180|; 
W. B. Kilgore, Phys. Rev. D 70 (2004) 031501 |hep-ph/0403128 



M. Grazzini and S. Frixione, hep-ph/0411399' 

[12] G. Heinrich, Nucl. Phys. Proc. Suppl. 116 (2003) 368 "hep-ph/0211144(; 

C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. D 69 (2004) 076010 |hep-ph/03ll3Tl| ; 
T. Binoth and G. Heinrich, Nucl. Phys. B 693 (2004) 134 hep-ph/04022l5"; 
G. Heinrich, Nucl. Phys. Proc. Suppl. 135 (2004) 290 hep-ph/0406332 . 

[13] A. Gehrmann-De Ridder, T. Gehrmann and G. Heinrich, Nucl. Phys. B 682 (2004) 265 
hep-ph/0311276 . 

[14] C. Anastasiou, K. Melnikov and F. Petriello, Phys. Rev. Lett. 93 (2004) 032002 hep-ph/0402280| . 
[15] C. Anastasiou, K. Melnikov and F. Petriello , prep-ph/0409Q'88| |hep-ph/050 1130 



[16] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Nucl. Phys. B 691 (2004) 195 
hep-ph/0403057 . 

[17] D.A. Kosower, Phys. Rev. D 57 (1998) 5410 hcp-ph/9710213 . 

[18] A. Gehrmann-De Ridder, T. Gehrmann and E.W.N. Glover, Nucl. Phys. Proc. Suppl. 135 (2004) 97 
|hep-ph/0407023| . 

[19] F.A. Berends and W.T. Giele, Nucl. Phys. B294 (1987) 700; 

M. Mangano, S. Parke and Z. Xu, Nucl. Phys. B298 (1988) 653; 
M. Mangano, Nucl. Phys. B309 (1988) 461; 

L.J. Dixon, Proceedings of "Theoretical Advanced Study Institute in Elementary Particle Physics (TASI 
'95): QCD and Beyond", ed. D. Soper, World Scientific (Singapore, 1996), p.539 |hep-ph/9601359| . 

[20] V. Del Duca, L.J. Dixon and F. Maltoni, Nucl. Phys. B 571 (2000) 51 hep-ph/9910563 ; 



F. Maltoni, K. Paul, T. Stelzer and S. Willenbrock, Phys. Rev. D 67 (2003) 014026 hep -ph/0209271| . 

[21] H.P. NiUes, Phys. Rept. 110 (1984) 1; 

H.E. Haber and G.L. Kane, Phys. Rept. 117 (1985) 75. 

[22] H.E. Haber and D. Wyler, Nucl. Phys. B 323, 267 (1989). 

[23] H. Bacr, X. Tata and J. Woodside, Phys. Rev. D 42 (1990) 1568. 



[24] J. Rosiek, Phys. Rev. D 41 (1990) 3464; |hep-ph/9511250 



12 



[25] S. Catani, Phys. Lett. B 427 (1998) 161 |hep-ph/9 802439l; 

G. Sterman and M.E. Tejeda-Yeomans, Phys. Lett. B 552 (2003) 48 |hep-ph/0210T30| . 



[26] A. Gchrmann-De Ridder, T. Gehrmann and E.W.N. Glover, hep-ph/0502110 



[27] F. Wilczek, Phys. Rev. Lett. 39 (1977) 1304; 

M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Phys. Lett. B 78 (1978) 443. 

[28] A. Gchrmann-De Ridder, T. Gehrmann and E.W.N. Glover, in preparation. 



13 



